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Abstract 



We present results from Experiment 864 for antiproton production and an- 
tideuteron limits in Au + Pb collisions at 11.5 GeV/c per nucleon. We have 
measured invariant multiplicities for antiprotons for rapidities 1.4 < y < 2.4 at 
low transverse momentum as a function of collision geometry. When compared 
with the results from Experiment 878 our measurements suggest a significant 
contribution to the measured antiproton yield from the decay of strange an- 
tibaryons. We have also searched for antideuterons and see no statistically 
significant signal. Thus, we set upper limits on the production at approxi- 
mately 3 x 10~ 7 per 10% highest multiplicity Au + Pb interaction. 
PACS numbers: 25.75.+r 
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I. INTRODUCTION 



Antimatter production in relativistic heavy ion collisions has been proposed as an excel- 
lent probe of the collision dynamics and possible phase transition to a quark-gluon plasma 
[IJ. The production of antiprotons at AGS energies (10-15 GeV/c) is near threshold in 
nucleon-nucleon collisions. Therefore multiple collisions, resonance excitations, and collec- 
tive effects may play a major role in significantly increasing the overall production rates [Q. 
Strange antibaryons, due to their larger mass and additional antistrange quark, are even fur- 
ther suppressed in initial nucleon-nucleon collisions. However, strangeness saturation and 
antimatter enhancement have long been predictions of quark-gluon plasma states. Thus, 
understanding the yields of non-strange and strange antibaryons is an important tool for 
distinguishing between various sources of enhanced production. 

Antibaryons have a large annihilation cross section (particularly at low relative mo- 
mentum), reaching levels of hundreds of millibarns. Thus, in these baryon rich colliding 
systems, there can be significant annihilation losses before they escape the collision region. 
The final experimentally measured yields represent the initial production minus the annihi- 
lation losses. The annihilation process in free space is well understood and experimentally 
parametrized; however, this annihilation may be modified in the dense particle environment 
where the initial attraction of baryon and antibaryon may be disturbed 0. 

Antideuteron production at AGS energies is actually below the energy threshold in single 
nucleon-nucleon collisions and thus antideuterons are expected to only be created through 
the coalescence of separately produced antiprotons and antineutrons which are close enough 
together in coordinate space and phase space at freeze-out. Because of the large energy 
required for their production, antideuteron yields are an excellent measure of the system's 
thermal temperature (assuming antimatter is equilibrated in these collisions). Also, coa- 
lescence yields would yield information about the spatial distribution of antinucleons (as a 
form of two particle correlation) jp. 

In the next section we describe the E864 spectrometer and the antiproton data sets. In 
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section three we present the measured antiproton invariant multiplicities and compare them 
with measurements made by other experiments. These comparisons lead us to consider the 
possibility of enhanced production of strange antibaryons in these collisions. Antimatter 
correlations in the form of events with two antiprotons and antideuteron production are 
discussed in section four. 

II. EXPERIMENT 864 
A. The E864 Spectrometer 

Experiment 864 was designed to search for novel forms of matter (particularly strange 
quark matter, or "strangelets" ) produced in heavy ion collisions at the Brookhaven AGS 
facility ||. In order to conduct this search, E864 has a large geometric acceptance and op- 
erates at a high data rate. A diagram of the spectrometer is shown in Figure Secondary 
particles produced from the Au + Pb reaction which are within the geometric acceptance 
traverse two dipole magnets (Ml and M2) and then multiple downstream tracking stations. 
Non-interacting beam particles and beam fragments near zero degrees pass above the down- 
stream spectrometer in a vacuum chamber, thus reducing interactions which would otherwise 
produce background hits in the detectors. The experiment does not measure at zero degrees 
(zero transverse momentum), but particles with at least 15 milliradians of angle pass through 
an exit window in the vacuum chamber. 

Charged particles are tracked using three time-of-flight (TOF) hodoscopes and two straw 
tube stations. The hodoscopes yield three independent charge measurements of dE/dx 
over the 1 cm thickness of the scintillator slats and provide three space-time points with 
time resolutions on the order of 120-150 ps. The straw stations provide more accurate 
position information for track projection back into the magnetic field region. Particles are 
identified by their time-of-flight from the target (yielding the velocity) and momentum. 
The momentum is determined by combining the charge measurement with the rigidity (R 
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= p/Z) from the track projection in the bend plane into the field region. The redundant 
measurements allow for excellent background rejection of ghost tracks and tracks originating 
from downstream interactions. 

A second particle identification measurement can be made using a hadronic calorimeter 
located at the end of the spectrometer || . The calorimeter consists of 754 individual towers 
constructed from lead sheets with scintillating fibers running almost parallel to the incident 
particle trajectory. The calorimeter yields good timing information (a ~ 400 ps for hadrons) 
and excellent hadronic energy resolution of 3.5% + 3A%/y/E (with E in GeV). For baryons, 
the calorimeter measures the particle's kinetic energy, which when combined with time- 
of-flight information gives a measure of the particle mass. For antibaryons, the energy 
measurement also includes the annihilation energy of the antibaryon and its annihilation 
partner. 

The experiment is able to perform high sensitivity searches by running at high rate with 
a special "late energy" trigger (LET) [[7]. The time and energy signals from each of 616 
fiducial calorimeter towers are digitized in flash ADC's and TDC's and used as inputs to a 
lookup table, which is programmed to select the particles of interest. Because there are many 
slow neutrons and many fast high energy protons, a simple time cut or an energy cut was 
determined to be insufficient for the trigger. The late energy trigger allows for the rejection of 
both of these abundant particles, while effectively triggering on slow (mid-rapidity) particles 
which deposit a large amount of energy. An antiproton of the same momentum as a proton 
or neutron will deposit an additional annihilation energy. If the Au+Pb interaction yields 
no towers firing the trigger, then a fast clear is sent out to the digitizers and the data is 
not recorded. The trigger yields an enhancement factor for antiprotons, antideuterons and 
strangelets of approximately 50 (under running conditions appropriate for each species). 

In order to determine the collision geometry (impact parameter) a charged particle mul- 
tiplicity counter is used. The E864 multiplicity counter |J is an annular piece of fast BC420 
scintillator placed around the beam pipe 13 cm downstream of the target and tipped at 
an angle of 8° to the vertical. It is 1 cm thick and subtends the angular range of 16.6° 
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to 45.0°. The annulus is separated into four quadrants and each quadrant is viewed by a 
photomultiplier tube. The total signal measured with this counter is proportional to the 
charged particle multiplicity of the collision. The integrated signal from the sum of the four 
quadrants is used to trigger on the centrality of the events by selecting events with a signal 
larger than a given threshold. 

B. Data Sets 

The data used in this analysis was collected in two separate data taking periods. Dur- 
ing the fall of 1995, the late energy trigger was strobed on the 10% most central Au+Pb 
interactions with the spectrometer magnets set for optimal acceptance for antideuterons 
(referred to as the "-0.75T" field setting). The LET curve was set to yield an enhancement 
factor of ~ 50 for antideuterons and negative strangelets. The data set includes over 90 
million recorded events, which effectively sampled approximately six billion central interac- 
tions. From this sample, over 50,000 antiprotons were identified. The mass distributions of 
antiprotons from a single rapidity and transverse momentum bin are shown in Figure 0. 

In the fall 1996 run, the LET was strobed on minimum-bias (93% of the geometric cross 
section) Au+Pb interactions and the LET curve and the spectrometer magnets were set 
for optimal antiproton acceptance (referred to as the "-0.45T" field setting). The LET 
yielded an enhancement factor for antiprotons > 50 under these conditions. However, in 
order to use the trigger effectively, the region of the calorimeter dominated by neutrons from 
the interaction had to be excluded. This reduced the geometric acceptance by roughly a 
factor of two. The data sample included 45 million recorded minimum bias interactions and 
approximately 50,000 antiprotons. These data samples represent the largest statistics for 
antiprotons produced in heavy ion collisions at the BNL-AGS. 

In both data sets, the beam momentum was measured using the E864 spectrometer 
magnets and a downstream beam counter located in the beam dump. The beam momentum 
of 11.5 GeV/c per nucleon was consistent with the beam momentum reported at extraction 
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from the accelerator once energy losses due to material in the E864 beam line were properly 
accounted for. The Au beam was incident on a 30% Pb target for the 1995 data set, while 
a 10% Pb target was used in 1996. 

III. ANTIPROTON INVARIANT MULTIPLICITIES 

A. E864 Measurements 

In E864 we explicitly measure the yield of antiprotons per Au+Pb interaction as a func- 
tion of centrality, and thus we directly calculate the invariant multiplicities. The invariant 
multiplicity for antiprotons is determined as follows: 

1 d 2 N _ 1 N counts 1 

27rp t dydpx 27rp7AyAp T N sampled e detec t x 

^accept X €trigger 

The total number of antiprotons N counts is determined in each separate bin in phase space 
and divided by the total number of sampled Au+Pb interactions N samp i e d. The counted 
antiprotons include only those antiprotons which fired the LET. Since the detector does 
not measure all the antiprotons produced in a given region of phase space, the invariant 
multiplicity must be corrected for the missed particles. These missed antiprotons are the 
result of the experiment's finite acceptance and various tracking and triggering efficiencies. 
The acceptance e accepi and detection efficiency edetect are calculated using a GEANT || 
simulation of the experiment in conjunction with real data. This simulation also included 
losses due to antiproton annihilation in the target as part of the acceptance correction. The 
production of antiprotons due to reinteraction of particles from the primary interaction with 
target nuclei was also considered and found to be negligible. 

The LET efficiency ^trigger is determined in each kinematic bin. This efficiency is deter- 
mined in one of two ways: for antiprotons measured in the 1995 ("-0.75T") 10% central data 
where the efficiency was somewhat low, a sample of antiprotons that did not fire the trigger 
was used to determine the efficiency. In the 1996 ("-0.45T") data, where the LET curve was 
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set for higher efficiency (~ 75%), the efficiency was determined from a Monte Carlo of the 
calorimeter response. 

The data from E864 is measured in a range of 50 < px < 275 MeV/c (where the limits 
are a function of rapidity) . The invariant multiplicities measured in E864 are approximately 
flat over the px range measured, as shown for the 1996 data in Figure [| Over such a 
small range in transverse momentum, the invariant multiplicities are not expected to change 
significantly. If, for example, the spectra follow a Boltzmann distribution, 

1 d 2 N 



2np T dydpT 



oc mre t b (2) 



(where = yp T + m 2 ), then with a temperature parameter of 200 MeV the invariant 
multiplicity at px = is only 6% higher than at px = 150 MeV/c. For comparison with 
other experiments, in each rapidity bin all the invariant multiplicities as a function of px 
are fit to a constant level. This level is assigned as the invariant multiplicity at px ~ 0, and 
an additional 6% systematic error is assigned due to the px = extrapolation. It should be 
noted that strong radial flow could affect the px = extrapolation as well. We feel that this 
effect should be within the estimated systematic uncertainty since the E864 data are quite 
flat as a function of px down to 50 MeV/c at midrapidity. 

The antiproton invariant multiplicities for 10% most central Au + Pb collisions from 
the 1995 data are given in Table | ||10|| . It should be noted that the statistical error in 



this data set is dominated by the contribution from the trigger efficiency (due to counting 
antiprotons which did not fire the trigger). The systematic error in the 1995 data (exclusive 
of the 6% due to the px = extrapolation) is estimated to be 15%, and is dominated by 
the uncertainty in the correction for the LET trigger efficiency. Systematic uncertainties 
also arise from our knowledge of the experimental acceptance (including the effect of the 
collimator in the first spectrometer magnet), track quality cuts, and the loss of tracks due 
to overlapping hits in the hodoscopes. 

For the 1996 data, the late-energy trigger was strobed on a minimum-bias sample of 
events selected by the multiplicity counter. The resulting multiplicity counter ADC dis- 
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tribution is shown in Figure f|. When selecting minimum bias events, it is important to 
consider the effect of interactions of the beam that do not occur in the target. Using special 
empty-target runs we have found that non-target interactions contribute less than 10% of 
the multiplicity distribution at low multiplicity, while the late-energy trigger further reduces 
this contamination to below 1% (see Figure |]). For the 1996 data, the antiproton invariant 
multiplicities are determined for different regions of the minimum bias multiplicity: 100-70% 
of the full distribution, 70-30%, 30-10% and 10%. These regions are shown in Figure ||. It 
is important to note that the LET rejection is a strong function of the multiplicity, and this 
must be properly accounted for when calculating the normalization in each centrality bin. 

The antiproton invariant multiplicities for the four multiplicity regions used in the 1996 
data set are listed in Table ||. In addition, the full minimum bias invariant multiplicity 
from the 1996 data set is also listed in Table |I[ The systematic error in these data points 
are estimated to be 10% (again, exclusive of the 6% previously described due to the px = 
extrapolation). As in the 1995 data, the systematic uncertainty in the 1996 data is also 
dominated by the uncertainty in the determination of the LET trigger efficiency. However, 
the size of the correction is smaller for the 1996 data due to the overall higher efficiency of 
the trigger setting. 

Figure P shows the 1995 ("-0.75T") and 1996 ("-0.45T") antiproton invariant multiplici- 
ties at pt = as a function of rapidity. The Gaussians shown are fits to the combined 1995 
and 1996 data, and are constrained to have a mean value at midrapidity (y = 1.6). There 
is excellent agreement between the 1995 and 1996 data in the rapidity range where the two 
data sets overlap. The rapidity widths measured in the data for the four multiplicty widths 
are 0.37±0.02 (100-70%), 0.41 ±0.02 (70-30%), 0.43±0.02 (30-10%) and 0.46±0.02 (10%), 
indicating a broadening of the rapidity spectrum at higher centrality. 

Also shown in Figure ^ (as open squares) are previously reported antiproton results from 
data taken in 1994 ("-0.45T") |TT|. It should be noted that the 1994 data is about 20% higher 
at midrapidity than indicated by the corresponding 1995 and 1996 data. This is within the 
statistical and systematic error previously quoted for the 1994 data. It is important to note 
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that the 1994 data was taken with an incomplete detector: two layers of S3 were missing 
along with a "plug" designed to reduce the occupancy in the downstream detectors due to 
interactions of beam-rapidity fragments with the vacuum chamber. The presence of the plug 
dramatically reduced the detector occupancy in the 1995 and 1996 data (and hence the size 
of the correction required for losses due to multiple hodoscope hits), and the presence of the 
additional S3 layers provided additional background rejection. 



B. Comparisons with Other Experiments 

Experiment 878 has measured antiproton yields as a function of collision geometry in 
reactions of Au+Au ions at 10.8 A GeV/c |12|,|13J. There are two differences between the 
reaction system between E878 and E864: (1) the target in E864 is Pb and (2) the beam 
momentum in E864 is higher at 11.5 A GeV/c. The target difference is quite small and is 
neglected in this comparison. However, the production of antiprotons is near threshold for 
nucleon-nucleon collisions at these energies and so the beam momentum difference must be 
accounted for. We assume that the ratio of antiproton yields in Au+Pb reactions at the two 
energies is proportional to the ratio of antiproton yields in p+p reactions at the two energies. 
Unfortunately, there is no usable data on antiproton production in p+p reactions covering 
this particular energy range. Therefore a parametrization of the production cross sections 
(derived from p+p data at higher energies and p+Be data from E802 at 14.6 GeV/c ||) is 
employed. Using this parameterization, one expects the ratio of the antiproton production 
cross sections at the two energies to be 1.5. The E878 invariant multiplicities are scaled 
up by this value. By considering fits to higher energy p+p data that do not include the 



E802 p+Be data at 14.6 GeV/c [1_4|, we estimate that this energy scaling contributes an 
additional 15% systematic error on the overall normalization of the scaled E878 points. 

Experiment 878 measures invariant multiplicities nominally at px = (which is really at 
transverse momenta less than ~ 30-50 MeV/c). Using the procedure previously outlined, 
we extrapolate the E864 measurements to pt = and compare the E864 and E878 mea- 
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surements in Figure |7|. While the two experiments agree well for low multiplicity collisions 
a substantial disagreement develops for more central collisions. For 10% central collisions, 
the E864 measurements at midrapidity are a factor of ~ 3.2 larger than the corresponding 
E878 data. 

It should be noted that both experiments do not use precisely the same definition of 
centrality: E864 measures the multiplicity of charged particles produced in the collisions, 
while E878 measures the 7 multiplicity (mostly from 7r° decay). In order to properly compare 
the two experiments the multiplicity ranges for both experiments must be converted to a 
(somewhat model-dependent) parameter. To do this, we have chosen to show the integrated 
antiproton yield at = versus the number of "first" nucleon-nucleon collions in each 
centrality range. In order to estimate the number of first collisions in each multiplicity 
range for the E864 data, a GEANT || simulation was used in conjunction with RQMD fj~8[ 



Au+Pb events to generate a trigger probability vs. impact parameter distribution for each 
multiplcicity region (see Figure |]). These distributions were then folded with distributions of 
the number of first collisions vs. impact parameter from a simple Glauber model calculation. 
A similar procedure was applied to the E878 data, using the results of a simulation of the 
E878 multiplicity array The results of this exercise (shown in Figure [II]) demonstrates 
that the E864 and E878 centrality ranges are quite similar. 

In Figure || we also compare measurements of the minimum-bias cross section for Au+Pb 



collisons at 11.5 A GeV/c with E878 and E886 fll6| . It should be noted that experiment E886 
only measured antiprotons from minimum bias collisions and thus there is no comparison 
as a function of centrality. As expected by the comparison of the E864 data with E878 for 
the four different centrality regions, the minimum bias invariant multiplicities measured in 
E864 are substantially larger than those measured by E878 and E886. 
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C. Strange Antibaryon Feed-Down 



There is a scenario which can reconcile the E864 and E878 results. Some of the antipro- 
tons measured by the various experiments may be the daughter product of weak decays 
of strange antibaryons (A, E+, E°, S°, S~, Q). This process is referred to as "feeding 
down" from the strange antibaryons into the antiprotons. Due to the significantly different 
designs of the two experiments, they have different acceptances from these decay product 
antiprotons. There are a number of antihyperon (Y) feed-down channels into the antiproton: 



A — > p+7T+ (65%B.F.) (3) 

TP -> A + 7^ p+7r + + 7 (100% x 65% B.F.) (4) 

E+^ p + 7i° (52% B.F.) (5) 

BP -> A + 7T ^ p+7r + + 7r (99% x 65% B.F.) (6) 

F -> A + 7T+^ p+7r + + 7r + (99% x 65% B.F.) (7) 



and multiple decay modes for the Q. The decay of the E° will produce additional A's which 
will be indistinguishable from those created in the primary collision. The decay of the A and 
the E + will produce p's whose production vertices do not coincide with the location of the 
primary interaction between the two nuclei. Therefore, the degree to which p's from these 
decays contribute to a measurement of p production will vary among experiments. 

Due to its large acceptance, the E864 spectrometer will detect p's from Y decay. E864 
does not have sufficient vertical resolution to distinguish p's from Y decay based on the 
vertical projection of a particle to the target, and the analysis cuts do not preferentially 
reject antiprotons from Y decay. Therefore, the p's detected in E864 are a combination of 
primary p's and p's from Y decay, in a ratio that reflects their production ratio. 
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The E878 collaboration have also evaluated the acceptance of their spectrometer for feed- 
down from Y decay. At midrapidity the acceptance for p's from A and S° decay is 14% of 
the spectrometer acceptance for primordial p's, and 10% of the p acceptance for E+ decays 
fT3|j ; the acceptance grows at higher rapidity. In what follows, we assume a uncertainty of 
±1% in the E878 acceptances for feed-down. 

Since both E878 and E864 measure a different combination of primordial p production 
and feed-down from Y decay, we can in principle separate the two components if we make 
two explicit assumptions: both E864 and E878 understand their systematic errors, and the 
entire difference between the two experiments can be attributed to antihyperon feed-down. 
It is important to note that in energy scaling the E878 results we have implicitly assumed 
that the K's scale with energy in the same way as the p's. 

Given an understanding of the errors involved, we can perform a statistical analysis of 
the Y/p ratio as a function of the E864 and E878 measurements (see Figures [11] and |I~2| ). 
This analysis results in the following limits on the ratio of Y/p : 



Y\ /A + £° + l.lE H 



P Jl'to V P 



0.02 (100-70%) 
0.10 (70-30%) 
> (98'X C,L.)\ 1.0 (30-10%) (8) 
2.3 (10%) 

0.2 (minimum bias) 

while the most probable value of this ratio is ~ 3.5 for 10% central collisions. The factor 
of 1.1 multiplying the £+ arises due to the different branching ratio and acceptance for 
the S+ compared to the A. The probability distributions in Figure |ll] were generated 
using the measured E864 and E878 invariant multiplicities for each centraility bin. The 
statistical errors on these measurements were treated as Gaussian, while systematic errors 
on the measurements, energy scaling, p? = extrapolation, and the E878 acceptance for 
feed-down were treated as defmining a limit around the measured values. 

E878 has not explicitly calculated their experimental acceptance for the doubly strange 
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H and the Q, and thus they are not explicitly included in the above formula. These heavier 
strange antibaryons are generally thought to be further suppressed and thus a small contri- 
bution. However, in light of the unexpectedly large "feed-down" contributions from strange 
antibaryons, one should be careful not to neglect their contribution to this ratio. 

This comparison indicates a Y/p ratio in Au+Pb collisions that is significantly greater 
than one at midrapidity and px = 0. It should be noted that if the F's and the p are produced 
with different distributions in y and pt, then the ratio of integrated yields of these particles 
will differ from the ratio at central rapidity and p? = 0. Preliminary results from Si+Au 
collisions based on direct measurements of p and A production by the E859 collaboration 
also indicate a ratio of integrated yields greater than one [pTSfl . For comparison, the A/p ratio 
in pp collisions at similar energies is ~ 0.2 |[20|| . 

An enhancement of antihyperons arises naturally in models that include a QGP, and 
therefore enhanced antimatter and strangeness production [[21], J] . Thermal models that 



use a temperature and baryon chemical potential derived from measured particle spectra 



also indicate that the Y/p ratio could be larger than one [22|. However, extremely large 
values of the Y jp are difficult to achieve in a thermal model unless the freezout temperature 
and/or K + /K~ are pushed beyond experimentally observed values. Transport models such 
as RQMD |TB[ predict the Y/p ratio to be less than one. Including in a cascade model 
conversion reactions such as 

p + K + -> vr + A (9) 

and a lower annihilation cross section for the A relative to the p enhances the Y/p ratio 
substantially p3f . However, such a model does not reproduce the trend with centrality seen 



in the E864 data. 



IV. ANTIMATTER CORRELATIONS 
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A. Double Antiproton Events 



In the large sample of events from the 1995 ( "-0.75T" ) run with a single antiproton within 
the experimental acceptance, there are some events with two identified antiprotons in the 
same event. These two antiproton events give insight into the possible correlated production 
of antimatter. Since the number of individual nucleon-nucleon collisions in each Au + Pb 
collision is large, if the sample of central events are similar in nature, the production of one 
antiproton should have very little relation to the production of a second antiproton. 

In the 1995 data set there are approximately 43,000 antiprotons with rapidity less than 
2.2, which were considered for this study. After corrections for background contributions, 
we find there are 3.8 events with two antiprotons. 

If we assume that the production of antiprotons is uncorrelated, we can calculate the 
number of two antiproton events expected. One can think of the nucleus-nucleus collision as 
many (n) nucleon-nucleon collisions each with a probability (p) of producing an antiproton. 
Since the probability (p) is small and the number of collisions (n) is large, we calculate the 
probability of producing two antiprotons in the same event using Poisson statistics. The 
probability of producing one antiproton is: 

Prob(l) = Ratesingies = n x p (10) 
The probability of producing two antiprotons is: 

Prob (2) = Rat6singlcs2 (11) 

Since we have measured the rate of single antiprotons into our detector Ratesingies, we calcu- 
late the expected number of two antiproton events at 1.8. The 90% confidence level upper 
limit on this number is five, which includes the experimentally measured value. 

Given the agreement with the assumption of uncorrelated production, there are limits 
we can set on the possible correlated production of antimatter. We postulate that there are 
two distinct classes of events within the 10% central Au+Pb sample: One class of purely 
hadronic reactions and one class with the formation of the quark-gluon plasma (QGP). 
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In Figure the predicted number of two antiproton events as a function of the fraction 
of QGP events Jqcp and the antimatter enhancement factor e is shown. The area in the dark 
box is where the predicted number of two antiproton events is greater than five and thus 
ruled out by the data at the 90% confidence level. If the QGP enhancement factor is small, 
the two antiproton yield is not changed significantly. Also, if most of the events are QGP, 
then regardless of the enhancement, there is no predicted increase in the two antiproton 
yield. However, if there is a large enhancement (e > 10) and the fraction of QGP events is 
between 5% to 25%, the yield of two antiprotons is significantly increased. These specific 
scenarios are ruled out by this measurement at the 90% confidence level. 

B. Antideuteron Search 

We have performed a search for antideuterons using the 1995 data set of central Au + 
Pb interactions taken at the "-0.75 T" magnetic field setting optimized for the acceptance 
of antideuterons. After processing the data, any tracked particle of charge negative one, 
rapidity y < 2.4, passing all track quality \ 2 selections, and having a reconstructed mass in 
the range 1.3 < m < 5 GeV/c 2 is considered a possible antinuclei candidate. The mass 
distribution of these candidates is shown in Figure |14]. The distribution is well fit by an 
exponential and has no significant signal at the antideuteron mass m = 1.874 GeV/c 2 . 

The experiment is able to reduce this background through an energy measurement us- 
ing our full coverage hadronic calorimeter. The calorimeter measures the deposited kinetic 
energy of hadrons in addition to the annihilation energy for antibaryons. The background 
processes expected to create high mass candidates in the tracking reconstruction are the 
result of neutrons which charge exchange in the vacuum exit window or air and produce a 
forward going proton traversing the downstream spectrometer. The protons have reasonable 
rapidity values, but reconstruct to erroneously large rigidities resulting from the incorrect 
assumption that the particle originated at the target. These candidates should leave signif- 
icantly less energy in the calorimeter than expected if they are actually protons compared 
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with real antideuterons or heavier antinuclei. 

Since these candidates are all assumed to be antimatter, the reconstructed calorimeter 
mass must account for the annihilation contribution. In studies of antiproton showers from 
test beam data and from the 1995 data, it was observed that only « 84% of the annihilation 
energy was recorded in the calorimeter. Thus, the calorimeter mass formula is modified to 
reflect this loss: 

maSS = ^68- (12) 

The tracking mass resolution is ~ 5%, which yields a a m = 0.094 GeV/c 2 for an- 
tideuterons. The distribution of calorimeter masses for candidates whose tracking mass is 
within ±2<j m of the antideuteron mass (1.687 < m < 2.061 GeV/c 2 ) is shown in Figure [15. 



One can see the peak mean value is less than 0.938 GeV/c 2 . Protons have a lower calorime- 
ter mass (< 0.938 GeV/c 2 ) when calculated using Equation [12| since they do not deposit 



any energy beyond their kinetic energy (there is no annihilation energy contribution). The 
background candidates appear to be protons as expected from charge exchange background. 
Most protons should not reconstruct such a large antimatter mass and fire the late-energy 
trigger. However, the calorimeter energy response has a non-Gaussian high side tail. These 
candidates are protons which occupy the high side tail part of the energy response distribu- 
tion. As can be seen in the plot, the calorimeter is a powerful tool for rejecting this proton 
background. A cut is then placed on calorimeter mass being greater than 1.600 GeV/c 2 . 

If one assumes that all of the observed candidates are from charge exchange background 
(really protons striking the calorimeter), then the background shape can be fit. The tracking 
mass distribution with no cut on the calorimeter mass is fit to a simple exponential func- 



tion as shown in Figure [T^. If the candidates are all protons striking the calorimeter, the 
calorimeter mass distribution should be the same regardless of the tracking mass. Thus, one 
can use the exponential function fit parameters from the tracking mass distribution with no 
calorimeter cuts to describe the tracking mass distribution with a calorimeter mass cut. 



The tracking mass distribution is plotted in Figure [LJ with a cut on the calorimeter 
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mass greater than 1.600 GeV/c 2 . There are ten candidates within the ±2cr m range of the 



antideuteron. The background fit distribution shown in Figure [H] is renormalized to the 
total number of counts and plotted. The exponential fit seems a reasonable description of the 
distribution. The total number of counts from the fit in the region of the antideuteron (within 
±2<r m ) is 9.0. Thus, there is no significant signal above background for the antideuteron. 

One can then ask, how many real antideuterons would there have to be to make a 
statistically significant peak above the background distribution. There are nine predicted 
background events, and thus the Poisson statistics 90% confidence level upper limit is 14.2. 
If more than 14.2 candidates were observed within the antideuteron mass range, there is 
less than a 10% chance that it is due to a statistical fluctuation in the background events. 
Thus, we set the 90% confidence level upper limit on antideuteron production at N Poisson — 

N Background 5.2. 

In order to translate this Poisson statistics limit into a total upper limit on the production 
of antideuterons, the various acceptances and efficiencies must be known. It is also possible 
using a specific production model to set the 90% confidence level upper limit on the invariant 
multiplicity in a given region of momentum space. In the discussion that follows, we will 
asuume a model in which the production in p T and rapidity (y) can be factored as: 

1 <lX - A Qe -2pT/<PT> e -(y-ycm) 2 /2o- y 2 ^ 13 ^ 



2irpT dydpx 

We have assumed a rapidity width a y = 0.5 and a mean transverse momentum < pr >= 1-00 
GeV/c. Using this production model, we set a 90% confidence level upper limit on the 
production of antideuterons at 2.78 x 10 -7 per 10% central Au + Pb interaction. 

One can relate the limit over all phase space to the limit on the invariant multiplicity 
(A ) at midrapidity and pt = 0. 

a _ -^TotalLimit / , . \ 

A ° = ( 27 r)3/V^ (M) 
The upper limit on the invariant multiplicity at midrapidity y = 1.6 and pt = is 1.4 x 
10~ 7 GeV~ 2 c 2 . We have tested the model dependency of these upper limits and find that 



with extreme ranges of production models, one can vary the upper limits by approximately 
±50%. 

Using our antiproton measurements and these upper limits, we calculate the 90% con- 
fidence level upper limits on the coalescence scale factor B 2 for antideuterons. This scale 
factor may be a function of where in momentum space the measurement is made, thus we 
give the limit at midrapidity (y = 1.6) and p t = 0. Our measured invariant multiplicity for 
antiprotons is 1.16 x 10 -2 GeV _2 c 2 (from the combined 10% central 1995 and 1996 data). 
The upper limit on the invariant multiplicity for antideuterons is 1.41 x 10~ 7 GeV~ 2 c 2 . The 
upper limit on the scale factor is: 

1 'i- v , j, 

< 1.0 x 10~ 3 GeVV 2 (15) 



Bo 



2-npt dydp. 



-(d) 



1 cPN 



2-npt dydpt 

This upper limit is shown as an arrow in Figure [T^, along with a comparison to coalescence 
scale factors measured at Bevelac and SPS energies p^|. 



This scale factor is significantly below the global value of 1.2 x 10~ 2 GeV 2 c -2 pre- 
dicted by the "simple" coalescence model. However, since this prescription has failed to 
describe systems where the collision volume is expected to be large compared with the 



deuteron/antideuteron size j27j, it is not surprising that it is in disagreement with the value 
obtained here. 

If the source distribution of antinucleons has a similar spatial extent as the nucleon source, 
then the scale factor for deuterons B 2 is expected to be the same as for antideuterons B 2 . 
Recently, E864 has presented measurements of protons and deuterons around midrapidity 
and low transverse momentum. The scale factor from the analysis of E864 light ion data 



is also shown in Figure 16. 



B 2 = 1.1 ± .4 x 10~ 3 GeVV 2 (16) 

The uncertainties are dominated by systematic errors in the deuteron and proton invariant 
multiplicities. This measured scale factor is at the same level as the upper limit for the 
antideuteron scale factor. We cannot determine whether B 2 is significantly lower than B 2 . 
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Thus, it is impossible to comment on whether the rate of antideuteron production is smaller 
due to preferential surface emission of antimatter. 

If we consider the most probable value of the ratio (Y)/p = 3.5 for 10% central collisions, 
the primordial antiproton multiplicity at midrapidity and p t = should be a factor of ~ 3.3 
lower than measured in E864. In this picture, the 90% confidence level upper limit on the 
coalescence scale factor would be 

B~ 2 < 1.1 x 10~ 2 GeVV 2 (17) 

This value is approximately at the level measured in p + A collisions, in accord with the 
"simple" coalescence model, where the collision volume is expected to be quite small. Thus, 
if the Y production is correctly calculated, the limit set on antideuteron production is not 
very significant in the context of coalescence models. 

There have been two previous measurements of antideuteron production in heavy ion 
collisions. The first was from the E858 experiment which observed two antideuterons in 
minimum bias Si + Au collisions at 14.6 A GeV/c [29|. They calculated a coalescence 



factor of approximately B 2 < 1.0 x 10~ 3 GeV 2 c~ 2 . While this value is consistent with 
our observation, it is difficult to make any direct comparison since the E858 value is for 
minimum bias collisions involving a much smaller projectile. The second measurement is 
from experiment NA52 in Pb + Pb central collisions at 160 A GeV/c |30| . They observe a 
coalescence factor of approximately B 2 ~ 5.0 ±3.0 x 10 -4 GeV 2 c -2 . They also find the factor 
for deuterons B 2 is the same within statistical and systematic errors. While our upper limit 
on antideuterons is consistent with their value, our deuteron coalescence factor is somewhat 
higher. This observation is not surprising due to larger source dimensions in the higher 
energy collisions studied by NA52. 



V. CONCLUSIONS 



We have presented results from Experiment 864 for antiproton production and an- 
tideuteron limits in Au+Pb collisions at 11.5 GeV/c per nucleon. We have measured in- 
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variant multiplicities for antiprotons above midrapidity and at low transverse momentum 
as a function of collision geometry. These measurements are within systematic errors of our 



previously reported results [11] , and, when compared with the results from Experiment 878, 
may indicate a significant contribution to the measured antiproton yield from the decay of 
strange antibaryons. 

We have also studied correlated production of antimatter using events with more than 
one antiproton and a search for antideuterons. For antideuterons we see no statistically 
significant signal. We set upper limits on the production at approximately 3 x 1CT 7 per 10% 
highest multiplicity Au+Pb interaction. 
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TABLES 



rapidity 


1 dN , 
2ix'Pt dydpr 


itp T = (xl0~ 2 , i 


n GeV 2 c 2 ) 


1.8 < y < 2.0 




1.00 ±0.05 




2.0 < y < 2.2 




0.67 ±0.03 




2.2 < y < 2.4 




0.36 ±0.02 





TABLE I. Antiproton invariant multiplicities at p T = from the 1995 (-0.75T) 10% central 
data. The errors listed are statistical only. 







1 dN at 

2vrp T dydp T 


PT = (xlO^ 2 , 


in GeV- 2 c 2 ) 




rapidity 


100-70% 


70-30% 


30-10% 


10% 


min. bias 


1.4 < y < 1.5 


0.12 ±0.02 


0.56 ±0.06 


0.87 ±0.08 


0.94 ±0.11 


0.48 ± 0.02 


1.5 < y < 1.6 


0.14 ±0.01 


0.55 ±0.02 


0.98 ± 0.04 


1.17 ±0.06 


0.51 ±0.01 


1.6 < y < 1.7 


0.13 ±0.01 


0.53 ±0.02 


0.90 ± 0.03 


1.15 ±0.05 


0.48 ±0.01 


1.7 < y < 1.8 


0.12 ±0.01 


0.52 ±0.02 


0.80 ±0.02 


1.11 ±0.03 


0.46 ±0.01 


1.8 < y < 1.9 


0.11 ±0.01 


0.44 ±0.02 


0.75 ± 0.02 


1.00 ±0.04 


0.39 ±0.01 


1.9 < y < 2.0 


0.093 ± 0.005 


0.38 ± 0.02 


0.68 ±0.02 


0.88 ±0.03 


0.34 ±0.01 


2.0 < y < 2.1 


0.060 ± 0.004 


0.30 ±0.01 


0.51 ±0.02 


0.71 ±0.03 


0.27 ±0.01 



TABLE II. Antiproton invariant multiplicities at pt = vs. centrality from the 1996 (-0.45T) 
minimum-bias data. The data are listed by multiplicity regions used in the analysis. The invariant 
multiplicities for the full minimum-bias sample are also listed. The errors listed are statistical only. 
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FIG. 1. The E864 detector in plan and elevation views, showing the dipole magnets Ml and 
M2, hodoscopes (HI, H2 and H3), straw tube arrays (S2 and S3) and hadronic calorimeter (CAL). 
The vacuum chamber is not shown in the plan view. 
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Mass (GeV/c 2 ) Mass (GeV/c 2 ) 

FIG. 2. Antiproton mass distributions for Z = — 1 particles from the 1995 ("-0.75T") data set. 
The left panel show the mass distribution for particles that triggered the LET, while the right panel 
shows the mass distribution for particles that did not. These mass distributions are for particles 
with rapidity 1.8 < y < 2.0 and transverse momentum 125 < pt < 150 MeV/c. 
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FIG. 3. Invariant multiplicities for antiprotons as measured in the four different centrality bins 
from the 1996 data. Note that the data are approximately flat in each rapidity interval over the 
measured range in transverse momentum. 
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FIG. 4. Multiplicity counter ADC distributions for INTO (minimum bias) and LET triggers for 
10% Pb and empty targets; the distributions are scaled to the same number of incident Au ions. 
For LET triggers the empty target contribution is less that 1% of the distribution for the lowest 
multiplicities. 
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FIG. 5. Multiplicity counter ADC distributions (corrected for empty target contributions) 
showing the centrality cuts used in the analysis. The effective LET rejection factor for each 
multiplicity region is also shown. 
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FIG. 6. Antiproton invariant multiplicities as a function of centrality showning the 1994, 1995 
and 1996 data sets from E864. The error bars shown are statistical only. Systematic errors are 
estimated to be 20% on the 1994 data, 15% for the 1995 data, and 10% for the 1996 data, not 
including a 6to pt = 0. The fits are constrained to have a mean value at midrapidity (y = 1.6). 
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FIG. 7. Antiproton invariant multiplicities at px = compared to E878 as a function of cen- 
trality. Note that while the two experiments agree for peripheral collisions, they disagree markedly 
for more central collisions. The E878 data have been scaled up as described in the text to account 
for the lower beam momentum. 
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FIG. 8. Antiproton invariant multiplicities from the 1996 (-0.45T) minimum bias data. Also 
shown is minimum bias data from E878 (scaled up to 11.5 GeV/c) and data from E886 (shown 
reflected about midrapidity, y = 1.6). 
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FIG. 9. Trigger probability vs. impact parameter for the E864 multiplicity counter in the four 
centrality regions used in the analysis of the 1996 data. These distributions were generated using 
RQMD events in a GEANT simulation of the E864 multiplicity counter. 
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FIG. 10. Integrated antiproton yield at px = vs. the mean number of first collisions. The 
number of first collisions are generated using simulations of the detector response and a Glauber 
model; the horizontal bar for each multiplicity range indicates the RMS of the first collisions 
distribution. 
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FIG. 11. Probability distributions for the antihyperon/antiproton ratio as a function of cen- 
trality. The distributions are generated from the ratio of the E864 antiproton measurements at 
midrapidity vs. the E878 measurements using the statistical and systematic errors. See the text 
for details. The 98% confidence levels are marked in for each distribution. 
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FIG. 12. 98% confidence level values for the antihyperon/antiproton as a function of centrality. 
The nominal pp level at 12 GeV/c is marked. Note the dramatic rise in the ratio as a function of 
collision centrality. 
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FIG. 13. The number of expected double antiproton events as a function of the fraction of total 
events which have a QGP formation and the relative antimatter enhancement factor in these QGP 
events to normal hadronic events. 
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FIG. 14. The upper panel shows the antideuteron candidate tracking mass distribution with 
no requirement on the calorimeter mass. The bottom panels shown the same distribution after a 
cut requiring the mass determined from the calorimeter cluster be larger than 1.6 GeV/c 2 . The 
vertical lines show the antideuteron mass window used in the analysis. 
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FIG. 15. Calorimeter reconstructed mass for antideuteron candidates with tracking mass 
1.687 < m < 2.061 GeV/c 2 . The line shows the mass cut at 1.6 GeV/c 2 applied in the anal- 
ysis. 



40 



^ 10 



o 



> 

CD 

o 

" -2 

m 1 o 



10 



-4 



Bevelac AA 



O O Q 



AGS AA 



SPS FNAL pA 



10 



-1 



E864 d/p 2 



v 



E864 90% U.L a/p : 



A=2 



10 



10 10^ 
T/Ap (GeV) 



FIG. 16. Coalescence scale factors as a function of kinetic energy per nucleon of the colliding 
beam. The upper limit on B2 from E864 is shown as an arrow. The E864 value for deuterons B2 
is also shown. The simple coalescence level is shown as a line with low energy AA and high energy 
pA results. 
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